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^ Abstract 
-i— > ■ 

In [Journal of Statistical Planning and Inference (141) (2011) 3697- 
3704], Roberto Fontana offers an algorithm for obtaining Sudoku ma- 
trices. Introduced by Geir Dahl concept disjoint pairs of S-permutation 
matrices [Linear Algebra and its Applications (430) (2009) 2457-2463] 
■ is used in this algorithm. Analyzing the works of G. Dahl and R. 

i Fontana, the question of finding a general formula for counting disjoint 

^sO ] pairs of n 2 x n 2 S-permutation matrices as a function of the integer n 

naturally arises. This is an interesting combinatorial problem that de- 
serves its consideration. The present work solves this problem. To do 
\ that, the graph theory techniques have been used. It has been shown 

that to count the number of disjoint pairs of n 2 x n 2 S-permutation 
matrices, it is sufficient to obtain some numerical characteristics of 
the set of all bipartite graphs of the type g = {R g U C g ,E g ), where 
j_j 1 V = R g U Cg is the set of vertices, and E g is the set of edges of the 

graphs, R g nC g = ®, \R g \ = \C, 



n. 
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1 Introduction 

Let Pij, 1 < i, j < n be n 2 square n x n matrices, which elements belong to 
the set 2) n 2 = {1,2,..., n 2 }. Then n 2 x n 2 matrix 
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is called Sudoku Matrix, if every row, every column and every submatrix Pjj, 
1 < hj < n make permutation of the elements of set 2) ra 2, i.e. every number 
s G {1,2, ... ,n 2 } is found just once in every row, every column and every 
submatrix Pjj. Submatrices Pij are called blocks of P. 

Sudoku is a very popular game. On the other hand, it is well known that 
Sudoku matrices are special cases of Latin squares in the class of gerechte 
designs [TJ. 

A matrix is called binary if all of its elements are equal to or 1. A 
square binary matrix is called permutation matrix, if in every row and every 
column there is just one 1. 

Let denote by E n 2 the set of all permutation n 2 x n 2 matrices of the 
following type 

" A u An ■■■ A ln 
A 2 \ A 2 2 ■ ■ ■ A 2n 



A 



A n \ A. 
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where for every s,t G {1,2,..., 72} A st is a square n x n binary submatrix 
(block) with only one element equal to 1. The elements of £ n 2 will be called 
S-permutation matrix. 

Two S-permutation matrices A and B will be called disjoint, if there are 



not i,j G 2) n 2 = {1,2, 
the condition = b i3 



, n } such that for the elements G A and 6y G B 
1 is satisfied. 



The following obvious proposition is given in [2]: 

Proposition 1 [2] A square n 2 x n 2 matrix P is Sudoku matrix if and only if 
there are n 2 mutually disjoint S-permutation matrices A ly A 2 , . . . , A n 2, such 
that P can be given in the following way: 

P = 1 • A 1 + 2 • A 2 H h n 2 ■ A n 2 



□ 
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In [5] Roberto Fontana offers an algorithm which randomly gets a family 
of n 2 x n 2 mutually disjoint S-permutation matrices, where n = 2,3. In n = 3 
he run the algorithm 1000 times and found 105 different families of nine 
mutually disjoint S-permutation matrices. Then according to proposition [1] 
he decided that there are at least 9! • 105 = 38 102 400 Sudoku matrices. 
This number is very small compared with the exact number of 9 x 9 Sudoku 
matrices. In [4J it has been shown that there are exactly 

9! • 72 2 ■ 2 7 ■ 27 704 267 971 = 6 670 903 752 021 072 936 960 

number of 9 x 9 Sudoku matrices. 

To evaluate the effectiveness of Fontana's algorithm is necessary to cal- 
culate the probability of two randomly generated matrices to be disjoint. As 
it is proved in [2j the number of S-permutation matrices is equal to (n\) 2n . 
Thus the question of enumerating all disjoint pairs of S-permutation matrices 
naturally arises. This work is devoted to this task. In the present paper, we 
will formulate and prove a formula (Theorem [1] and in another form Theorem 
[2]) for counting disjoint pairs of S-permutation matrices. To do that, we will 
use graph theory techniques. 

2 A formula for counting disjoint pairs of S- 
permutation matrices— the graph theory ap- 
proach 

Let us denote by S m the symmetric group of order m, i.e. the group of 
all one-to-one mapping of the set 3) m = {1,2,..., m} in itself. If x G 2) m , 
p G S m , then the image of the element x in the mapping p we denote by p(x). 
Bipartite graph is the ordered triplet 

g= (R g UC g ,E g ), 

where R g and C g are non-empty sets such that R g nC g = 0, the elements of 
which will be called vertices. E g C R g x C g = {(r,c) \ r G R g ,c G C g } - the 
set of edges. Repeated edges are not allowed in our considerations. 

Let g' = (R g > U C g >,E g r) and g" = (R g » U C g »,E g »). We will say that the 
graphs g' and g" are isomorphic and we will write g' = g" , if Rg' — Rg" > 
Cg> = C g " and there are p G S m and a G S n , where m = \R g >\ = \Rg"\ and 
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n = \C g '\ = \C g "\, such that (r, c) G E g < (p(r),a(c)) G -Eg». The object 
of this work is bipartite graphs considered to within isomorphism. 

For more details on graph theory see for example [3] or [6]. 

Let n and k be two positive integers and let < k < n 2 . Let us denote 
by <2»n,fc the set of all bipartite graphs without multiple edges of the type g = 
(R g U Cg,E g ), considered to within isomorphism, such that \R g \ = \C g \ = n 
and \E g \ = k. 

Let g = (R g U C g , E g ) G & n ,k f° r some natural numbers n and k and let 
v G V g = Rg U Cg. By 7 (v) we denote the set of all vertices of V g , incident 
with v, i.e. u G j(v) if and only if there is an edge in E g connecting u and 
v. If v is an isolated vertex (i.e. there is no edge, incident with v), then by 
definition ^(v) = and |7(i>)| = 0. If v G R g , then obviously 7(1;) C C g , and 
if v G C g , then 7(f) C R g . Obviously 

J2h(v)\ = 2k. 

VGVg 

Let g = (R g U C g , E g ) G nj fc and let u, v G V g = R g U C s . We will say 
that u and v are equivalent and we will write u ~ v if 7(14) = 7(1?). If it 
and f are isolated, then by definition it ~ v if and only if u and d belong 
simultaneously to R g , or C s . The above introduced relation is obviously an 
equivalence relation. 

By Vgj^ we denote the obtained factor-set (the set of the equivalence 
classes) according to relation ~ and let 

V ah = {A 1 ,A 2 ,...,A S }, 

where A« C R g , or A* C C g , i — 1, 2, . . . s, 2 < s < 2n. We put 

Si = I Af|, 1 < 5i < n, i—l,2,...,s 

and for every g G we define multi-set (set with repetition) 

[9] = {Si, 6 2 , ---Ss}, 

where 5i,5 2 , ■ ■ ■ ,S S are natural numbers, obtained by the above described 
way. Obviously 

s 

i=l 
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\iz\Z% ... z n is a permutation of the elements of the set 2) n = {1,2, ... ,n} 
and we shortly denote p this permutation, then in this case we denote by p{i) 
the i-th element of this permutation, i.e. p(i) — Zi, i — 1,2, ... ,n. 

Lemma 1 The number of all n x n binary matrices with exactly k, k = 
0,1, ... ,n 2 ones, denoted by b(n, k) is equal to 



Proof. Let A = [ciij] nxn be nxn binary matrix with exactly k ones. Then 
we construct graph g = (R g U C g , E g ), such that the set R g = {r\, r 2 , . . . , r n } 
corresponds to the rows of A, and C g = {ci, c 2 , . . . , c n } corresponds to the 
columns of A, however there is an edge connecting the vertices r$ and Cj 
if and only if = 1. The graph, which has been constructed, obviously 
belongs to (Sy^. 

Conversely, let g = (R g U C g , E g ) e & n ,k- We number at a random way 
the vertices of R g with the help of the natural numbers from 1 to n without 
repeating any of the numbers. We analogously number the vertices of C g . 
This can be made by (n!) 2 ways. Then we construct binary nxn matrix 
A = [dij] , such that = 1 if and only if there is an edge in E g connecting 
the vertex with number i of R g with the vertex with number j of C g . As 
9 £ ®n,k, then the matrix, that has been constructed, has exactly k units. It 
is easy to see that when q,r 6 2) ra , g-th and r-th rows of A are equal to each 
other (i.e. the matrix A does not change if we changes the places of these 
two rows) if and only if the vertices of R g corresponding to numbers q and r 
are equivalent according to relation ~. Analogous assertion is true about the 
columns of the matrix A and the edges of the set C g , which proves formula 



Corollary 1 For every positive integer n, n > 2 and every natural number 
k such that < k < n 2 the following formula is true: 





se[g] 



□ 





□ 
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Let II n denotes the set of all n x n matrices, constructed such that tt G II n 
if and only if the following three conditions are true: 

i) the elements of n are ordered pairs of numbers where 1 < i, j < n; 

ii) if 

b x ) (02,62) •■■ (a n ,b n )] 

is the z-th row of n for any % G 2)„ = {1,2, ...,n}, then the numbers 
a\ a 2 ... a„ in this order construct permutation of the elements of the 
set 2)„ 

iii) if 

(aiA) 
(02,62) 

(«n, 6 n ) 

is the j-th column of ix for any j G 2) n , then numbers 61, 62, . . . , 6 n in this 
order construct permutation of the elements of the set 2) n . 

From the definition it follows that every row and every column of any 
matrix of the set IT n can be identified with permutation of elements of the 
set 2) n . Conversely for every (2n)-tuple ((pi, p 2 , . . . , p„), (ci, <r 2 , . . . , <r n )), 
where p; = pi(l) p»(2) . . . Pi(n), aj = aj(l) aj(2) . . . (Tj(n), 1 < i,j < n are 
permutations of elements of 2)n, then the matrix 



7T 



( Pi (1),cti(1)) (Pi(2),a 2 (l)) 
(p 2 (l),a 1 (2)) (p 2 (2),a 2 (2)> 



(p n (l),o"i(n)) (p„(2),a 2 (n)) 
is matrix of IL . Hence 



(pi(n),a n (l)) 
(p 2 (n),cr n (2)) 

(p n (n),tf n (n)) 



|n n | = (n! 



,2n 



(5) 



We say that matrices tt', 7t" G IT n , where 7r' = [»'•„•] , 7r" = \p"^\ 
are disjoint, if p'^ ^ p" ^ for every pair of indices i,j G 2) n . 

Lemma 2 Le£ n be a positive integer, n > 2. T/ien t/iere one to one 
correspondence between the sets E n 2 and LT n . 



Proof. Let A G S n 2. Then A is constructed with the help of formula (J2J) 
and for every z,j G 2) n in the block there is only one 1 and let this 1 has 
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coordinates (cii,bj). For every i,j G 2)„ we obtain ordered pairs of numbers 
(a*, 6j) corresponding to these coordinates. As in every row and every column 
of A there is only one 1, then the matrix [c%] nx „) where = (ai,bj), 
1 < hj < n i which is obtained by the ordered pairs of numbers is matrix of 
Il n , i.e. matrix for which the conditions i), ii) and iii) are true. 

Conversely, let [a i:j ] nxn G Il n , where a y = (a^bj), i,j G 2) n , a*, bj G 2) n . 
Then for every i,j G 2) n we construct binary n x n matrices with only 
one 1 with coordinates (a,i,bj). Then we obtain the matrix of type: 
According to the properties i), ii) and iii), it is obvious that the obtained 
matrix is S-permutation matrix. 

□ 

From Lemma H] and formula (jSJ) it follows the proof of Proposition [TJ in 

0- 

Corollary 2 (Proposition [T] of [2J) The number of S-permutation matrices 
of order n 2 is given by 

|S„ 2 | = N) 2n (6) 

□ 

It is easy to see that with respect of the described in Lemma [2] one to one 
correspondence, every pair of disjoint matrices of E n 2 will correspond to a 
pair of disjoint matrices of Il n and conversely every pair of disjoint matrices 
of Il n will correspond to a pair of disjoint matrices of S n 2. 

Corollary 3 The number of all pairs of disjoint matrices of S n 2 is equal to 
the number of all pairs of disjoint matrices ofU n . 

□ 

Lemma 3 The number q(n,k) of all ordered pairs of matrices (n',7r"), where 
7r', 7t" G Il n and having at least k, k = 0, 1, . . . , n 2 equal elements is equal to 

II (n-l7(«)l)! 
— - — ~Yj ( ) 

Proof. Let n' = Tp-.l , ix" = \p''A G Yi n and let tt' and 7r" have 
exactly k equal elements (by components). Then we definitely obtain the 



q(n,k) = (d) 2(n+1) V 
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binary n x n matrix A = [a-ij} nxn , such that = 1 if and only if p\- = 
p'-j, i,j G 2) n . According to Lemma [T] there is only one graph g G & n ,k 
corresponding to the matrix A and respectively to the ordered pair of matrices 
(tt', tt") G n x n. 

Conversely, let g = (R g UC g , E g ) G and let V g = R g UC g . Then there 



are — — — in number n x n binary matrices, corresponding to g, following 

the described in Lemma [1] rule and let A = [aij] nxn be one of them. Let 
7r = [Pij] nxn be an arbitrary matrix of n n . We will find the number h(jr, A) 
of all matrices tt' = [pL] nXn G Il n , such that p'^ = Pij if = 1. It is 
allowed to exist s, t G 2) n such that a st = and p' st = p st . Let the i-th. row 
(z = 1, 2, . . . , n) of 7r correspond to the permutation of the elements of 
2)„ and let i-th row of the matrix A correspond to the vertex G i? 9 of the 
graph g (see Lemma d]). Then there are (n— |7(Vi)|)! in number permutations 
of elements of 2) n , such that if a# = 1, then pj(t) = p-(t), t G 2)„. 
Analogously we can prove the corresponding assertion for the columns of 
7r. Hence, h(ir, A) = j J ( // — |t( v ) I)'- From all the things that have been 

described above, it follows that for every tt G Il n there are 

matrices of Il n , which have at least k elements, which are equal to the cor- 
responding elements of re (by components). But according to formula ([5]) 
|n n | = (T2,!) 2n , from where it follows formula ([7]). 

□ 

The purpose of the present work is to prove the following theorem: 

Theorem 1 Let n > 2 be a positive integer. Then the number D n 2 of all 
possible ordered pairs of disjoint matrices in S n 2 is equal to 

n (»-i7(*)D! 

D n2 = (n\r + (n\)^ f^i-l)" VeR9U ° 9 • (8) 

Se[g] 
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The number d n 2 of all non- ordered pairs of disjoint matrices in £ n 2 is 
equal to 

1 

d n 2 = -D n 2 (9) 

Proof. Let n be a positive integer, n > 2. Then applying Lemma [2] 
and the inclusion-exclusion principle we obtain, that the number D n 2 of all 
possible ordered pairs of disjoint matrices of £ n 2 is equal to 

n 2 

D n 2 = |n n | 2 + ]T(-l) fc 4(n,fc), 

k=l 

where the function q(n, k) is defined in Lemma [3] and may be calculated with 
the help of formula ([7]), and |n„| with the help of formula ([5]). Then we 
obtain the proof of formula ([8]) . 

Having in mind that the "disjoint" relation is symmetric and irreflexive 
we obtain formula (Q. 

Theorem [T] is proved. 

□ 

Let n and k are positive integers and let g G <S n ,k- We examine the 
ordered {n + l)-tuple 

*(g) = (Mg),Mg),---iM9)), (10) 

where ipi(g), i = 0,1, ... ,n is equal to the number of vertices of g incidents 

n 

with exactly i number of edges. It is obvious that ii^i(g) = 2A; is true for 

i=l 

all g G & n ,k- Then formula ([8]) can be presented 



D n 2 = {n\f n + (n!) 2 (™ +1 ) ^ - 



n-iW i{9) 



k=i g e<s n>k II ^ 
se[g] 



Since {n — n)\ — 0! = 1 and [n — (n — 1)]! = 1! = 1, then 



n-2 



n [(» 

D n2 = (n\r + (n\) 2(n+1) ^(-l) fc £ t= ° • (11) 

fc=l S6®„,fe H ^ ! 

«6fa] 



Consequently, to apply formula ( ITT]) for each bipartite graph g e <& n <k and 
for the set (S n ^ of bipartite graphs, it is necessary to obtain the following 
numerical characteristics: 



n-2 



<5e[ 9 ] 



and 



;i2) 



(13) 



9(n,k) = J2 <%) 

Example. When n = 3 and k = 6 the set 3j6 = {g x , g 2 , g 3 , g<t, gs, 9&\ 
consists of six bipartite graphs, which are shown in Figure [TJ 









Figure 1: The set of bipartite graphs 
For graph g\ 6 3j6 we have: 

[ 9l ] ={1,1,1,1,1,1} 



w(gi 



*(^i) = (0,0,6,0) 

[(3-0)!]° [(3-1)!]° 
1! 1! 1! 1! 1! 1! 
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For graphs g 2 G (83,6 and g 3 G <5 3 £ we have: 

M = [^3] = {1,1,1,1,2} 



¥(02) = *(<&) = (0,1,4,1) 

[(3-0)!]°[(3-l)!] 1 60-2 1 
= wOft.) = 1!2 , = -g- = 1 

For graph p 4 G 63,6 we have: 

[^4] = {1,1, 1,1, 1,1} 

^(^ 4 ) = (0,2,2,2) 

[(3-0)!]°[(3-l)!] 2 6°-2 2 
^ " 1!1! 1! 1! 1! 1! " ~T ~ 4 
For graphs #5 G 03,6 and g§ G (83,6 we have: 

= [g 6 ] = {1,2,3} 



*(^ 5 ) = *(^ 6 ) = (1,0,3,2) 

[(3-0)!] 1 [(3-l)!]° e 1 ^ 1 
"(flO = "foO = ^ = W = 2 

Then for the set <5 3i6 we get: 

0(3,6)= ^ u(g) = 1+1+1+4+^+^=8 

9e<»3,6 

□ 

All bipartite graphs with n = 2 and n = 3 together with the above nu- 
merical characteristics are described in [8]. 

Using the numerical characteristics (|T2|) and (|T3|) . we obtain the following 
variety of Theorem [TJ 



11 



Theorem 2 

n 2 

D n , = (n!) 4 " + (n\) 2 ^ £(-l)*0(n, fe), (14) 

fe=i 

where 9(n, k) is described using formulas fTB) and fTB). 

□ 

3 Conclusions and future work 

In order to apply Theorem [1] or Theorem [2] it is necessary to describe all 
bipartite graphs to within isomorphism g = (R g U C g ,E g ), where \R g \ = 
\C g \ = n. For n = 2 and n = 3 the application of Theorem [2] is described in 
[S] in full details. According to calculations, we obtain: 

• The number _D 4 of all disjoint ordered pairs of S-permutation matrices 
in n = 2 is equal to 

D 4 = 144 (15) 

• The number of all disjoint non-ordered pairs of S-permutation ma- 
trices in n = 2 is equal to 

d 4 = \d± = 72 (16) 

• The number D$ of all disjoint ordered pairs of S-permutation matrices 
in n = 3 is equal to 

A) = 1 260 085 248 (17) 

• The number d 4 of all disjoint non-ordered pairs of S-permutation ma- 
trices in n = 3 is equal to 

d 9 = ^D 9 = 630 042 624 (18) 

Problem 1 Let n > 2 is a natural number and let G be a simple graph 
having (n\) 2n vertices. Let each vertex of G be identified with an element of 
the set S n 2 of all n 2 x n 2 S-permutation matrices. Two vertices are connected 
by an edge if and only if the corresponding matrices are disjoint. The problem 
is to find the number of all complete subgraphs of G having n 2 vertices: 
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Note that the number of edges in graph G is equal to d n 2 and can be 
calculated using formula (JSJ and formula Q (respectively formulas (II 2p . 
(USD, JUD and (J9D). 

Denote by z n the solution of the Problem [1] and let a n is the number of all 
n 2 x n 2 Sudoku matrices. Then according to Proposition [1] and the method 
of construction of the graph G, it follows that the next equality is valid: 

(19) 



(n 2 ) 



We do not know a general formula for finding the number of all n 2 x n 2 
Sudoku matrices for each natural number n > 2 and we consider that this 
is an open combinatorial problem. Only some special cases are known. For 
example in n = 2 it is known that o~ 2 = 288 [7j. Then according to formula 
(USD we get: 

a 2 288 
4! 24 

In [4] it has been shown that in n = 3 there are exactly 

a 3 = 6 670 903 752 021 072 936 960 = 

= 9! x 72 2 x 2 7 x 27 704 267 971 = 
2 20 x 3 8 x 5 1 x 7 1 x 27 704 267 971 1 ~ 6.671 x 10 21 
number of Sudoku matrices. Then according to formula (fT9j) we get: 



03 = 6 670 903 752 021 072 936 960 = 
9! 362 880 
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